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A  classical  nova  occurs  when  material  accreting  onto  the  surface  of 
a  white  dwarf  in  a  close  binary  system  ignites  in  a  thermonuclear 
runaway1,2.  Complex  structures  observed  in  the  ejecta  at  late  stages3"5 
could  result  from  interactions  with  the  companion  during  the  common- 
envelope  phase6,7.  Alternatively,  the  explosion  could  be  intrinsically 
bipolar,  resulting  from  a  localized  ignition  on  the  surface  of  the  white 
dwarf  or  as  a  consequence  of  rotational  distortion9,10.  Studying  die 
structure  of  novae  during  the  earliest  phases  is  challenging  because 
of  the  high  spatial  resolution  needed  to  measure  their  small  sizes1 '. 
Here  we  report  near-infrared  interferometric  measurements  of  the 
angular  size  of  Nova  Delphini  20 1 3,  starting  one  day  after  the  explo¬ 
sion  and  continuing  with  extensive  time  coverage  during  the  first 
43  days.  Changes  in  the  apparen  t  expansion  rate  can  be  explained  by 
an  explosion  model  consisting  of  an  optically  thick  core  surrounded 
by  a  diffuse  envelope.  The  optical  depth  of  the  ejected  material  changes 
as  it  expands.  We  detect  an  ellipticity  in  the  light  distribution,  sug¬ 
gesting  a  prolate  or  bipolar  structure  that  develops  as  early  as  the 
second  day.  Combining  the  angular  expansion  rate  with  radial  veloc¬ 
ity  measurements,  we  derive  a  geometric  distance  to  the  nova  of 
4.54  ±  0.59  kiloparsecs  from  the  Sun. 

Nova  Delphini  2013  (also  known  as  V339  Del)  was  discovered12  by 
K.  Itagaki  on  14  August  2013  Universal  Time  (ut)  at  14:01.  We  began 
an  intensive  observing  campaign  to  measure  the  size  of  the  nova  with  the 
CHARA  Array13,  an  optical/infrared  interferometer  located  on  Mount 
Wilson,  California.  Our  observations  began  within  1 5  h  of  the  discov¬ 
ery  and  within  24  h  of  the  detonation  itself.  We  measured  the  size  of 
the  expanding  ejecta  as  the  nova  rose  to  peak  brightness  and  continued 
monitoring  it  for  a  total  of  27  nights  between  UT  2013  August  15  and 
September  26. 

We  measured  the  angular  diameter  of  Nova  Del  2013  by  filtinga  uni¬ 
formly  bright  circle  to  the  visibility  amplitudes  of  the  interference  fringes 
recorded  during  each  night  (Extended  Data  Table  1 ).  We  plot  the  expan¬ 
sion  curve  in  Fig.  1 .  The  dotted  line  shows  a  linear  fit  to  the  angular  diam¬ 
eters  during  the  first  27  days  after  the  explosion.  The  inset  panel  shows 
an  apparent  deceleration  during  the  first  four  nights  compared  with  lines 
of  constant  velocity.  During  the  last  week,  lire  measurements  show  a 
large  jump  in  the  effective  size  of  the  nova  compared  with  an  extrapo¬ 
lation  of  the  linear  fit.  For  comparison,  we  show  the  visible  and  infrared 
light  curves  in  Fig.  2. 

The  apparent  changes  in  slope  of  the  expansion  curve  can  be  explained 
by  a  two-component  model  consisting  of  an  optically  thick  pseudo¬ 
photosphere  surrounded  by  an  optically  thin  halo.  We  approximated 


the  intensity  distribution  projected  on  the  sky  using  a  uniformly  bright, 
circular  core  surrounded  by  a  circular  ring.  We  fitted  the  two-component 
model  to  the  interferometric  visibilities  and  minimized  the  total  y~  simul¬ 
taneously  across  all  nights.  We  fixed  the  time  of  detonation  at  modified 
Julian  day  t0  =  M]D  56518.277,  computed  by  extrapolating  the  first 
two  pre-discovery  photometric  measurements14  back  to  the  quiescent 
flux1546.  From  our  interferometric  data,  we  measured  an  expansion  rate 
of  0.156  ±  0.003  milli-arcseconds  (mas)  per  day  for  the  core  diameter 
and  a  size  ratio  of  1 .73  ±  0.02  between  the  outer  ring  and  core.  We  allowed 
the  flux  ratio  between  the  ring  and  core  to  vary'  on  a  nightly  basis. 

Figure  3  shows  how  the  percentage  of  light  from  the  ring  changes  over 
time.  During  the  first  two  nights,  the  core  and  ring  have  similar  surface 


Figure  1  |  Expansion  cum  of  Nova  Del  2013.  The  angular  diameters  were 
measured  by  filling  a  circular  disk  to  the  interferometric  data.  The  dotted  line 
shows  a  linear  fit  for  days  0-27.  The  inset  panel  zooms  in  during  the  first  week 
and  shows  dotted  lines  with  velocities17  of  613kms-i  and  2,500  km  s-1  at  a 
distance  of  4.54  kpc.  The  apparent  deceleration  during  the  first  week,  along  with 
tire  jump  in  size  during  the  last  week,  can  be  explained  by  a  two-component 
model  consisting  of  a  circular  core  surrounded  by  a  ring  where  the  flux  ratio 
changes  over  time.  The  grey  region  shows  the  expansion  rate  of  the  core 
flower  edge)  and  the  ring  (upper  edge).  The  solid  line  shows  the  effective  size 
obtained  by  fitting  the  two-component  model  visibilities  as  a  single  circular 
disk.  Error  bars  represent  la  uncertainties  computed  from  a  bootstrap  analysis. 
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Figure  2  |  Infrared  light  curee  of  Nova  Del  2013.  The  magnitudes  and  error 
bars  were  computed  from  the  mean  and  standard  deviation  of  counts  recorded 
on  the  detector  during  each  interferometric  observation  collected  at  the 
CH  ARA  Array.  The  data  were  obtained  in  the  H-band  using  two  different  beam 
combiners  (blue  squares  and  green  diamonds)  and  in  the  K’-band  (red 
triangles).  For  comparison,  we  plot  daily  averages  of  photometric 
measurements  in  the  V  band  downloaded  from  the  American  Association  of 
Variable  Star  Observers  (AAVSO;  black  circles).  The  rise  to  peak  brightness  was 
slower  in  the  infrared  compared  to  the  visible. 

brightnesses.  Therefore,  at  the  earliest  times,  the  nova  can  be  approxi¬ 
mated  by  a  single  uniform  disk  component  where  the  optically  thick 
pseudo-photosphere  extends  to  the  outermost,  fastest-moving  layers 
of  the  ejecta.  After  the  peak  in  the  visible  light  curve,  the  amount  of  flux 
in  the  ring  drops,  which  indicates  that  the  outer  layers  became  optic¬ 
ally  thin  and  the  pseudo-photosphere  moved  towards  the  inner,  slower- 
moving,  and  denser  layers.  Therefore,  the  apparent  deceleration  in  the 
expansion  curve  during  the  first  few  nights  is  probably  caused  by  a  dimin¬ 
ishing  contribution  from  the  outer  layers.  Over  the  next  few  weeks,  tire 
absolute  flux  of  both  components  decreased;  however,  the  flux  from 
the  core  dropped  at  a  steeper  rate  compared  with  the  ring.  Therefore, 
the  percentage  of  light  in  the  ring  increased  relative  to  the  core.  These 
changes  are  consistent  with  the  spectral  evolution'7  of  Nova  Del  2013 


Time  since  outburst  (days) 


Figure  3  |  Changes  in  the  flux  ratio  of  the  two-component  model.  The  solid 

line  shows  a  fit  where  the  ratio  changes  smoothly  over  time.  For  the  first  three 
nights,  the  ring  represents  the  outer  boundary  of  the  optically  thick  pseudo¬ 
photosphere.  After  the  peak  in  the  visible  light  curve,  the  flux  in  the  ring  drops. 
On  nights  4-27,  the  percentage  of  light  from  the  ring  increases  at  an 
approximately  linear  rate.  On  the  last  four  nights,  the  ring  contributes  an 
average  of  68%  of  the  total  light.  Error  bars  represent  lcr  uncertainties  derived 
from  a  least-squares  lit  to  the  visibility'  data. 


(see  http://www.astrosurf.com/aras/novae/Nova2013Del.html).  Optical 
spectra  initially  showed  P  Cygni  absorption  features  during  the  optically 
thick  fireball  stage,  which  disappeared  after  about  five  days.  Afterwards, 
strong  emission  lines  began  to  dominate  the  spectrum,  which  are  thought 
to  form  in  the  outer,  optically  thin  layers2.  The  rising  strength  of  the  emis¬ 
sion  lines  relative  to  the  continuum  was  also  seen  in  spectra  at  infrared 
wavelengths  {Extended  Data  Fig.  1 ).  In  the  two-component  model,  we 
assumed  a  constant  ratio  between  the  size  of  the  core  and  the  halo.  In 
reality',  changes  in  the  optical  depth  of  the  ejecta  are  likely  to  be  more 
complex  and  could  be  investigated  further  by  fitting  a  physical  model18,19 
simultaneously  to  interferometric  and  spectroscopic  data.  Understand¬ 
ing  how  the  optical  depth  changes  within  the  expanding  ejecta  affects 
the  physical  interpretation  of  novae  light  curves20. 

During  the  last  five  nights,  the  two-component  model  shows  an 
enhancement  in  the  flux  of  the  outer  layers,  with  the  ring  contributing 
about  68%  of  the  total  light.  This  coincides  with  an  increase  in  lire  infrared 
flux  from  the  nova21  (see  Fig.  2).  This  suggests  the  formation  of  dust  in 
the  outer  layers,  perhaps  in  denser  clumps  within  the  ejecta2,  that  con¬ 
tributes  thermal  emission.  Alternatively,  as  the  ejecta  expands,  a  recom¬ 
bination  front  might  propagate  through  the  material22.  Near-infrared 
free-free  emission  from  ionized  gas  in  the  inner  regions  would  decline 
and  shift  the  effective  size  boundary  to  the  outer  layers.  This  interpre¬ 
tation  is  consistent  with  the  larger  size  of  the  effective  photosphere  at 
later  dates,  but  it  does  not  explain  the  brightening  in  the  near-infrared 
light  curve.  The  effective  size  of  our  two-component  model  is  shown  as 
the  solid  line  in  Fig.  1 ;  the  model  reproduces  the  changes  in  slope  of  the 
expansion  curve. 

We  reconstructed  model- independent  images  of  Nova  Del  2013  using 
the  Markov  chain  imager25  (MACIM)  for  data  obtained  on  three  nights 
during  lire  first  week  with  sufficient  baseline  coverage.  As  shown  in  Fig.  4, 
the  images  showa  striking  similarity  to  tine  two-component  model,  with 
an  optically  thick  core  surrounded  by  a  halo  of  fainter  emission. 

Combining  the  angular  expansion  rate  with  the  radial  velocity'  of  the 
ejected  material  provides  a  way  of  measuring  the  distance2  to  Nova  Del  2013. 
The  velocilies  reported  in  Are  literature17  range  from  —600  kins'-1  to 
—  2,500  km  s-!,  consistent  with  the  expectation  that  lines  form  at  dif¬ 
ferent  layers  in  the  expanding  atmosphere  with  a  radial  outflow  velo¬ 
city  proportional  to  the  distance  from  the  white  dwarf  explosion  site. 
We  selected  the  Si  n  wavelength  k  =  6,347  A  and  2.  =  6,371  A  absorp¬ 
tion  lines  observed  during  the  first  week  to  represent  the  deeper  layers, 
because  they  probably  form  in  the  dense,  cooler  gas  immediately  above 
the  continuum-forming  layer.  From  an  analysis  of  spectra  downloaded 


Figure  4  |  Model  and  reconstructed  images  of  Nova  Del  2013.  The  two- 
component  model  consists  of  a  circular  core  surrounded  by  a  fainter  ring  (top  row). 
Tire  images  reconstructed  using  MACIM  show  a  similar  spatial  structure  (bottom 
row).  The  data  used  in  the  reconstructions  were  obtained  with  the  Michigan 
Infrared  Combiner  (MIRC)  at  the  CHARA  Array  on  UT  20 1 3  August  17, 19  and 
21  ((  =  3.0  days,  5.0 days  and  6.9  days  after  the  outburst)  and  have  good  baseline 
coverage  on  the  sky.  We  scaled  the  flux  by  the  infrared  H-band  magnitude 
measured  on  each  night  to  show  how  the  surface  brightness  changes. 
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from  the  archive  of  the  Astronomical  Ring  for  Access  to  Spectroscopy' 7, 
we  estimated  the  outflow  speed  near  the  continuum-forming  layer  to 
be  Vejertjon  =  613  ±  79kms~1 2 3 4 5 6.  Combining  this  velocity  with  the  angu¬ 
lar  expansion  rate,  we  derived  a  distance  of  4.54  ±  0.59  kpc  to  Nova 
Del  20 1 3.  Empirical  evidence  suggests  that  novae  have  prolate  spheroid24 
or  bipolar19  shapes,  so  the  true  distance  will  depend  on  its  geometry  and 
inclination  in  the  plane  of  the  sky25.  Assuming  a  probable  inclination  of 
45c  for  Nova  Del  2013  (S.  Shore,  personal  communication,  16  April  2014), 
we  estimate  that  these  effects  contribute  to  a  systematic  difference  in 
the  distance  of  less  than  5%.  Using  the  light  curve  decay  times  and  the 
reddening  reported15  for  Nova  Del  2013  (assuming  an  Rv  =  3.1  extinc¬ 
tion  law),  live  maximum  magnitude  to  rate  of  decline  (MMRD)  relation26 
gives  a  distance  ranging  from  3.3  kpc  to  4.1  kpc,  depending  on  the  adop¬ 
tion  ofthe  linear  or  nonlinear  MMRD  formulations.  The  expansion  par¬ 
allax  also  confirms  the  distance  of  4.2  kpc  based  on  a  spectral  comparison27 
to  the  nova  OS  Andromedae  1986. 

To  investigate  the  development  of  asymmetric  features  in  Nova 
Del  2013,  we  fitted  a  uniformly  bright  ellipse  to  the  visibilities  on  nights 
with  coverage  over  a  range  of  position  angles  on  the  sky  (from  two  to 
nine  days  after  the  outburst;  Extended  Data  Table  2).  For  each  night,  the 
ellipse  produces  a  significantly  lower  compared  with  the  circular  disk 
(P<  0.01  for  most  nights;  see  Methods).  The  nova  appears  to  be  around 
1 3%  larger  along  the  major  axis  compared  with  the  minor  axis,  suggest¬ 
ing  that  the  ejecta  could  be  bipolar.  For  the  remaining  nights  with  limited 
sky  coverage,  we  fixed  the  axis  ratio  and  position  angle  to  their  average 
values  (dmajor/^minor =  1.13  ±  0.07,  PA  =128'  ±  30")  and  solved  forthe 
size  of  the  ellipse.  On  average,  we  find  that  the  circular  disk  diameters 
underpredict  the  mean  angular  diameter  of  the  ellipse  by  about  2%. 
Therefore,  the  assumption  of  circular  symmetry  in  the  two-component 
model  should  not  strongly  bias  our  results. 

In  addition  to  the  visibility  amplitudes  that  provide  information  on 
the  size  and  shape  of  the  source,  interferometric  closure  phases  obtained 
when  combining  the  light  from  three  or  more  telescopes  indicate  devia¬ 
tions  from  point  symmetry  (whether  the  brightness  is  the  same  when 
reflected  through  a  point  at  the  centre  of  the  source  distribution).  The 
closure  phases  from  the  first  two  sets  of  data  obtained  with  all  six  tele¬ 
scopes  (t  =  3.0  and  5.0  days)  are  consistent  with  zero  within  the  uncer¬ 
tainties,  suggesting  that  the  source  was  point-symmetric  within  the 
resolution  limits.  By  the  third  set  of  observations  with  all  six  telescopes 
(f  =  6.9  days;  see  Fig.  4),  the  closure  phases  rise  steadily  to  about  60' 
(Extended  Data  Fig.  2),  suggesting  the  detection  of  a  point  asymmetry 
in  the  brightness  distribution.  If  the  outflow  is  bipolar,  as  suggested  by 
our  elliptical  fits,  then  this  could  indicate  a  difference  in  the  brightness 
between  the  two  lobes  (perhaps  due  to  viewing  angle)  or,  alternatively, 
the  development  of  clumpy  structures  within  the  expanding  material. 

The  ellipticity  and  closure  phase  asymmetries  indicate  that  non- 
spherical  structures  developed  as  early  as  a  few  days  after  the  outburst. 
Along  with  results  from  previous  interferometric  studies  of  novae28-30, 
this  suggests  that  novae  explosions  might  be  inherently  bipolar  or  that 
the  elliptical  shape  develops  early  during  the  common  envelope  phase. 

Online  Content  Methods,  along  with  any  additional  Extended  Data  display  items 
and  Source  Data,  are  available  in  the  online  version  ofthe  paper,  references  unique 
to  these  sections  appear  only  in  the  online  paper. 
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METHODS 

Interferometric  observations  with  the  CHARA  Array.  The  CHARA  Array13  is 
an  optical/infrared  interferometer  located  on  Mount  Wilson,  California.  The  array 
has  six  1  m  telescopes  arranged  in  a  Y  configuration  with  baselines  ranging  from 
34m  to  331  m.  The  longest  baseline  provides  a  spatial  resolution  of  0.5 mas  at  a 
wavelength  of  1.6  )im.  We  observed  Nova  Del  20  i  3  in  the  infrared  part  of  the  spec¬ 
trum  using  the  CLASSIC  two-beam  combiner31,  the  CLIMB  three-beam  combiner31, 
and  the  MIRC  six  beam  combiner32.  The  data  were  reduced  using  the  standard 
reduction  pipeline  for  each  instrument.  Data  products  include  the  squared,  nor¬ 
malized  visibility  amplitude  of  the  interference  fringes  for  the  two-baseline  com¬ 
biners  and  the  squared,  normalized  visibility  amplitudes  and  closure  phases  for  the 
instruments  that  combine  the  light  from  three  or  more  telescopes  (see  Extended  Data 
Figs  2-4) .  Before  and  after  each  observation,  we  observed  unresolved  stars  or  stars 
with  known  angular  diameters33  35  to  calibrate  the  interferometric  measurements. 
Extended  Data  Table  1  lists  the  UT  date,  modified  Julian  Day  (MJD),  time  since  the 
outburst  (f0  =  MJD  565 1 8.277),  beam  combiner,  telescope  configuration,  filter  and 
calibrators  used  during  each  observation.  The  calibrated  data  will  be  available  through 
the  Optical  Interferometry  Database  developed  by  the  Jean-Marie  Mariotti  Center 
(http://www.jmmc.fr/oidb.htm).  On  five  nights  we  also  collected  data  using  two 
visible  light  beam  combiners;  these  data  will  be  discussed  in  a  subsequent  publi¬ 
cation  (O.C.,  D.M.,  P.T.,  V.M.,  I.T.-B.,  D.P.K.B.,  E.  Lagadec,  M.I.,  N.N.,  A.M,  P.S., 
R.T.Z.,  G.H.S,  T.t.B.,  D.R.G.,  H.A.M.,  T.B.,  C.D.F.,  N.V.,  N.S.,  J.S.  and  L.S.,  manu¬ 
script  in  preparation). 

In  addition  to  the  interferometric  measurements,  we  also  computed  the  infrared 
magnitudes  of  Nova  Del  2013  using  the  number  of  counts  recorded  on  die  detec¬ 
tor.  We  compared  the  counts  from  Nova  Del  2013  with  the  calibrator  stars  observed 
immediately  before  and  after.  We  calibrated  the  photometry  based  on  the  2MASS 
magnitudes36  of  the  calibrators.  We  present  the  photometry  in  Extended  Data 
Table  3.  Figure  2  shows  die  infrared  light  curve  of  Nova  Del  2013.  For  comparison, 
we  also  plot  daily  averages  of  V -  band  photometric  measurements  downloaded 
from  AAVSO5’.  The  rise  to  peak  brightness  was  slower  in  the  infrared  compared 
to  the  visible  part  of  the  spectrum.  Based  on  a  second-order  polynomial  fit  to  the 
Fl-band  and  K'-band  measurements  near  the  maximum  ( <4.3  mag),  the  near- 
infrared  light  curve  reached  a  maximum  of  3.61  mag  at  MJD  56524.6  ( t  =  6.3  days). 
The  maximum  brightness  in  the  visible15  occurred  at  an  earlier  time  of  MJD  56520,4 
(t  =  2.1  days)  at  V  =  4.46  mag. 

Correcting  for  the  effects  of  line  emission.  The  emergence  of  strong  emission 
lines  in  the  infrared  spectrum  of  Nova  Del  2013  (such  as  Br,'  in  the  K-band)  will 
broaden  the  central  envelope  of  the  interference  fringe  and  create  sidelobes  farther 
out51’.  W  e  investigated  the  effects  of  the  emission  lines  on  the  calibration  of  the  inter¬ 
ferometric  data  obtained  with  the  CLASSIC  and  CLIMB  beam  combiners.  The  reduc¬ 
tion  code  outputs  tile  mean  power  spectrum  for  each  target  We  fitted  a  rectangular 
top  hat  function  to  each  power  spectrum  and  calculated  the  mean  centre  and  width 
for  the  observations  on  each  night  By  comparing  the  results  from  the  nova  with  the 
calibrators,  we  found  that  the  central  wavelength  of  the  filter  remained  constant  over 
time.  We  also  found  that  the  effective  bandwidth  of  the  filter  remained  the  same 
during  the  first  two  weeks  and  then  decreased  by  about  14%  ±  3%  as  the  emission 
lines  appeared  in  the  nova  spectra  i  see  Extended  Data  Table  4).  Essentially,  the  emis¬ 
sion  lines  put  a  larger  fraction  of  the  flux  at  the  centre  of  the  wavelength  region, 
thereby  effectively  decreasing  the  full-width  at  half-maximum  (FWHM)  of  tire  filter. 
The  CLASSIC/CL1MB  reduction  software  assumes  a  constant  top-hat  function  for 
the  width  of  the  filter.  Therefore,  when  the  effective  bandwidth  decreases,  the  visi¬ 
bility  is  overestimated  and  the  angular  size  of  the  nova  is  underestimated. 

To  correct  for  this  effect,  we  multiplied  the  CLASSIC  CLIMB  visibilities  by  the 
ratio  of  the  mean  effective  width  of  the  nova  observations  relative  to  the  calibrators 
measured  on  each  night.  We  added  the  uncertainly  in  the  bandwidth  ratio  in  quad¬ 
rature  with  the  uncertainties  in  the  risibilities  reported  by  the  reduction  code.  We 
applied  this  correction  starting  on  UT  2013  August  28,  when  the  effects  of  the  emis¬ 
sion  lines  became  measurable.  During  the  first  two  weeks,  when  tlae  contribution 
from  the  emission  lines  had  a  negligible  effect  on  the  calibration,  we  assumed  the 
default  bandwidth  used  in  the  reduction  code,  but  included  the  average  uncertainty 
in  the  nova-to-calibrator  bandwidth  ratio  of  ±0.164  in  the  error  budget  so  that  all 
of  the  CLASSIC  and  CLIMB  visibilities  would  have  similar  weights  during  the 
global  two-component  fit.  For  the  uniform  disk  fits  to  the  data  from  each  night,  tire 
bandwidth  correction  led  to  the  measurement  of  the  angular  size  of  the  nova  being 
larger  at  the  approximately  5%  level  after  the  emission  lines  developed 
Uniform  disk  and  uniform  ellipse  fits.  For  each  individual  night  of  CHARA  data, 
we  measured  the  angular  diameter  of  Nova  Del  2013  by  fitting  a  uniformly  bright, 
circular  disk  to  live  visibilities  and  performing  a  Levenberg-Marquardt  least-squares 
minimization  using  the  IDL  mpfit  package36  (http://purl.com/net/mpfit).  The  uni 
form  disk  diameters  (0UD),  reduced  /;  (where  v  is  the  number  of  degrees  of  free¬ 
dom),  and  number  of  visibility  measurements  Ar{W)  on  each  night  are  listed  in  the 
last  three  columns  of  Extended  Data  Table  1 . 


For  five  nights  with  sufficient  baseline  coverage  over  a  range  of  position  angles, 
we  fitted  a  uniformly  bright  ellipse  to  the  visibility  data.  The  fitted  parameters  include 
the  size  of  the  major  and  minor  axes  and  position  angle  of  the  major  axis  (f)major, 
(9mlluln  PA);  the  values  are  given  in  Extended  Data  Table  2.  For  four  of  the  five 
nights,  performing  an  F-test  indicates  that  the  uniform  ellipse  produces  a  signifi¬ 
cant  improvement  in  >f  compared  with  a  uniform  circle,  at  a  significance  level  of 
P  =  0.01.  On  UT2013  August  16,  the  improvement  is  only  at  the  0.10  significance 
level. 

For  both  the  uniform  disk  and  uniform  ellipse  fits,  we  determined  uncertainties 
in  the  parameters  using  a  bootstrap  approach*.  In  each  iteration,  we  created  a 
synthetic  sample  of  risibility  measurements  from  our  original  set  of  measurements 
by  randomly  selecting  the  same  total  number  of  risibilities,  with  replacement.  There¬ 
fore,  in  the  synthetic  sample,  some  measurements  are  repeated  and  others  are  left 
out  W e  applied  normally  distributed  uncertainties  to  the  synthetic  measurements 
(using  llte  measured  uncertainties  from  theoriginal  data).  We  then  fitted  a  uniform 
circle  or  uniform  ellipse  to  tire  synthetic  data  set  to  determine  tire  best-fitting  param¬ 
eters  for  that  iteration.  We  performed  1 ,000  iterations  and  computed  uncertainties 
based  on  the  standard  deviation  ofthe  parameter  distributions.  These  uncertainties 
are  reported  in  Extended  Data  Tables  1  and  2. 

We  expected  the  size  of  the  nova  to  vary  over  the  course  of  the  night.  The  maxi¬ 
mum  length  of  observing  time  was  about  6  h  on  UT  201 3  August  1 6  and  September  6. 
Given  an  expansion  rate  of  0. 156  mas  per  day,  we  would  expect  the  size  of  the  nova 
to  change  by  around  0.04  mas  over  the  course  ofthe  observations.  For  the  three  nights 
with  lire  smallest  uncertainties  on  the  angular  diameter  measurements  (±0.005  mas 
on  UT  2013  August  1 7, 1 9  and  2 1 ),  the  length  of  observing  time  ranged  from  30  min 
to  45  min,  so  the  change  in  size  is  expected  to  be  <0.005  mas.  For  all  other  nights, 
the  length  of  observing  time  varied  from  about  30  min  to  3  h,  corresponding  to  a 
change  in  size  of 0.003  mas  to  0.02  mas.  For  these  nights,  the  quoted  uncertainties 
on  the  angular  diameters  in  Extended  Data  Table  1  are  larger  titan  ±0.02  mas,  which 
should  account  for  any  changes  in  size  over  the  course  of  the  night.  For  the  two- 
component  model  described  in  the  next  section,  we  fitted  a  dynamic  model  dit  ectiy 
to  the  visibilities  versus  time,  so  it  accounts  for  changes  in  size  over  the  course  of 
the  nighL 

Two-component  model.  The  two-component  model  consists  of  a  uniformly 
bright,  circular  disk  surrounded  by  a  uniformly  bright  ring.  The  angular  diameter 
of  the  central  uniform  disk  0UD  and  the  outer  diameter  of  the  ring  0rim,  are  given  by: 

0UD(O  =  §('~to)  (1) 

$ring(0  =  Qmg0UD(O  (“) 

dS 

where  —  is  the  angular  expansion  rate  of  the  central  core  and  Cr5ng  is  the  ratio  of  the 
outer  ring  diameter  relative  to  the  uniform  disk  diameter. 

We  determined  the  two  global  parameters  {—  and  Crjng)  by  minimizing  the 
cumulative  y~  between  the  measured  and  model  visibilities  across  all  nights.  To  do 
this,  we  searched  through  a  grid  of  values  for  the  expansion  rate  and  size  ratio  and 
calculated  the  size  of  each  component  during  the  times  of  observation.  Then,  for  each 
night  of  observation,  we  performed  a  Levenberg-Marquardt  least  squares  mini¬ 
mization  using  the  IDL  mpfit  package39  to  solve  for  the  best  fitting  flux  ratio  between 
die  ring  and  die  core  component. 

There  are  degeneracies  between  the  overall  size  of  the  nova  and  the  flux  ratio  of 
the  two  components;  a  smaller  expansion  rale  can  be  compensated  by  placing  more 
flux  into  the  outer  component.  Therefore,  we  followed  a  two-step  procedure  in 
fitting  the  expanding  two-component  model  simultaneously  to  all  nights.  Initially, 
we  allowed  the  flux  in  the  ring  relative  to  the  core  to  vary  each  night,  but  with  a  limit 
set  so  that  the  surface  brightness  of'  the  ring  would  not  exceed  that  of  the  core.  The 
initial  restriction  prevented  a  large  increase  in  the  fitted  brightness  of  the  ring  at 
later  nights.  Using  die  grid- search  described  above,  we  determined  an  expansion  rate 
of  0.156  ±  0.003  mas  per  day  for  the  core  diameter  and  a  size  ratio  of  1.73  ±  0.02 
between  the  outer  ring  diameter  and  die  core  diameter.  We  then  fixed  the  expan¬ 
sion  rate  and  size  ratio  and  solved  for  the  ring-to  core  flux  ratio  for  each  night  after 
removing  die  restriction  on  the  surface  brightness  on  all  but  the  first  three  nights 
(we  would  not  expect  a  clearing  of  the  central  region  at  such  early  times}.  This  allows 
for  a  moderate  brightening  of  the  ring  at  later  times. 

In  Fig.  1,  the  solid  line  shows  the  effective  angular  size  of  the  two-component 
model.  This  curve  was  generated  by  assuming  that  the  percentage  of  light  from  the 
ring  varies  smoothly  over  time  (solid  line  in  Fig.  3}  and  computing  the  model  visi¬ 
bilities  for  each  night  based  on  the  expansion  rate  and  size  ratio.  We  then  fitted  a 
single,  uniform,  circular  disk  to  the  model  visibilities  from  each  night  to  compute 
the  effective  size.  Extended  Data  Fig.  5  shows  how  the  size  and  flux  ratios  ofthe  two 
components  change  during  die  nights  of  observation. 
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Emission  line  (lux  from  infrared  spectroscopy.  To  investigate  how  the  emission 
line  flux  changed  in  the  infrared,  we  used  low  dispersion  spectra  obtained  with  the 
TripleSPEC  spectrograph41 42  on  the  200-inch  Hale  Telescope  at  Palomar  Obser¬ 
vatory  on  UT  2013  August  20  and  23.  We  reduced  the  spectra  using  the  Spexlool 
reduction  program  modified  for  TripleSpec43,  which  includes  routines  to  correct 
for  detector  artefacts44.  We  calibrated  and  removed  the  telluric  lines  in  the  nova 
spectrum  using  the  xtdlcor  program44,  which  compares  a  high  signal-to-noise  spec¬ 
trum  of  Vega  to  observations  of  a  nearby  AO  star,  SAO  88500,  taken  immediately 
after  each  observation  of  the  nova.  We  also  used  infrared  spectra46  obtained  at  Mount 
Abu,  India42,  on  UT  2013  August  28  and  29,  and  September  8  and  20. 

We  collected  these  spectra  into  two  regions  covering  the  wavelength  range  of  the 
CHARA  CLASSIC/CLIMB  H-band  and  K'  -band  filters  and  normalized  the  flux  to 
unity  at  1 .65  pm  and  2.2  pm,  respectively.  The  H-band  is  dominated  by  the  upper- 
level  transitions  of  the  hydrogen  Brackett  series  plus  a  number  of  C 1  lines,  and  the 
K'-hand  is  dominated  by  Br; ,  ITei  /  =  2.0585  pm  plus  some  Cl  lines47. 

To  estimate  the  flux  from  the  continuum,  we  made  spline  fits  to  the  lower  enve¬ 
lope  of  each  spectrum.  We  then  summed  both  the  continuum  flux  and  the  emis¬ 
sion  flux  (total  minus  continuum)  across  the  wavelength  bands  to  form  a  ratio  of 
the  integrated  emission  line  flux  relative  to  the  integrated  continuum  flux.  As  shown 
in  Extended  Data  Fig.  1,  the  emission  strength  grew  over  the  time  frame  of  the 
CHARA  observations  ofNova  Del  20)  3  and  began  to  decrease  in  the  final  spectrum. 
These  changes  are  similar  to  the  changes  we  measured  in  the  effective  bandwidth  of 
the  interferometric  observations. 

The  strengthening  in  the  emission  lines  helps  to  explain  why  the  expansion  curve 
ofNova  Del  2013  became  steeper  alter  the  first  week.  The  emission  lines  form  over 
a  greater  radial  extent  in  the  expanding  envelope  than  does  the  continuum  from 
the  pseudo  photosphere.  Thus,  if  the  outflow  follows  a  Hubble  type  relation  with 
the  velocity  proportional  to  radial  distance48,  then  the  emission  lines  sample  parts 
of  the  envelope  that  were  expanding  faster  than  the  photosphere.  We  interpolated 
linearly  between  the  spectroscopic  measurements  to  obtain  the  ratio  of  the  emis¬ 
sion  line  flux  to  the  continuum  flux  at  the  times  of  the  CHARA  observations  and 
used  this  as  the  initial  estimate  of  the  flux  ratio  in  the  halo  ring  relative  to  the  uni¬ 
form  disk  core  in  the  two-component  fit. 

The  drop  in  the  emission-to-continuum  flux  ratio  after  the  fourth  week  can  be 
allributed  to  the  onset  of  dust  formation.  This  can  be  seen  in  the  K-band  spectra, 
where  the  continuum  slopes  do™  in  the  first  spectra  and  becomes  flat  in  the  last 
spectrum.  This  is  consistent  with  flux  from  a  dust)'  component  that  contributes  more 
to  the  continuum  and  increases  with  wavelength.  Our  two-component  model  based 
on  the  interferometric  observations  suggests  that  if  dust  formation  occurred,  that  it 
happened  over  a  larger  spatial  scale  than  that  of  the  expanding  pseudo  photosphere. 

We  also  made  measurements  of  the  flux-weighted  central  wavelength  of  each 
infrared  spectr  um.  These  were  essentially  constant  to  within  a  per  cent  or  so  during 
these  observations,  consistent  with  our  results  from  the  analysis  of  the  fringe  power 
spectra. 

Measuring  the  radial  velocity  of  the  optically  thick  core.  To  compare  the  angu¬ 
lar  expansion  rate  with  the  radial  velocity  of  the  nova  outflow,  we  needed  to  decide 
which  features  in  the  optical  spectra  are  most  representative  of  the  kinematics  of 
the  near-infrared  continuum.  The  velocities  reported  in  the  literature17  range  from 
—600  km  s  1  to  —2,500  km  s~  consistent  with  the  expectation  that  lines  form  at 
different  layers  in  the  expanding  atmosphere  with  a  radial  outflow’  velocity  propor¬ 
tional  to  distance  from  the  white  dwarf  explosion  site.  The  optical  depth  unity  bound¬ 
ary  for  the  near  infrared  continuum  must  form  at  higher  densities  deep  inside  the 
outflow  (compared  to  the  P  Cygni  lines  and  emission  lines  that  form  further  out  in 
lover  density  gas).  The  absorption  lines  observed  during  the  first  week  are  prob¬ 
ably  the  most  representative  of  the  deeper  layers,  because  they  form  in  the  dense, 
cooler  gas  immediately  above  the  continuum  forming  layer.  We  selected  the  Si  II 
a  =  6,347  A  and  6,371  A  absorption  lines  for  measurement  because  they  were 
observed  throughout  the  first  week  and  are  relatively  free  of  blending  and  inter¬ 
ference  from  Earth’s  atmospheric  lines  and  nova  emission  lines. 

We  downloaded  six  spectra  with  a  high  resolving  power  (R  =  10,000)  obtained 
over  the  first  week  by  O.  Garde  from  die  archive  of  the  Astronomical  Ring  for  Access 


to  Spectroscopy17  and  transformed  them  to  a  unit  continuum  on  a  heliocentric 
wavelength  grid.  We  measured  radial  velocities  for  the  Si  n  lines  by  cross-correlation 
with  a  model  spectrum  (T^g  —  8,000  K,  logg  =  2.0,  solar  melallicity)  based  upon 
ATLAS  atmospheres  from  R.  Kurucz.  The  mean  velocity  and  standard  deviation 
from  this  sample  (MJD  56519.8  —56524.8)  are  —598  ±  51  km  s-1.  The  local  stan¬ 
dard  of  rest  radial  velocity'  at  the  probable  distance  of  the  nova  is  + 15  km  s  so 
the  implied  outflow  velocity  is  —613  ±  51  kms~\  However,  we  do  not  know  the 
peculiar  velocity  of  the  nova  relative  to  its  local  standard  of  rest,  so  we  applied  a 
representative  velocity  dispersion  of  60  km  s~ 1  for  white  dwarfs  to  revise  the  uncer¬ 
tainty  in  our  final  estimate  of  the  outflow  speed  near  the  continuum  forming  layer, 
Vtjection  =  613  ±  79 km  s~]. 

We  note  that  this  is  strictly  an  upper  limit  for  tire  velocity  of  the  continuum 
forming  layers  because  the  absorption  lines  form  above  the  continuum  in  a  some 
what  faster  outflow.  Additionally,  the  estimate  reflects  the  geometry'  of  the  outflow, 
which  is  dependent  on  the  axial  inclination  tor  a  bipolar  ejection.  For  example,  if 
the  absorption  lines  formed  in  a  spherically  expanding  photosphere,  then  we  might 
expect  to  observe  line  boundaries  ranging  from  zero  (flux  from  the  limb)  to  —  Voutflow 
(flux  from  the  centre),  with  a  mean  velocity  between  these  limits.  However,  the 
observed  Si  It  lines  are  rather  narrow  (FWHM  is  280  km  s -1),  contrary  to  this  pre¬ 
diction.  On  tile  other  hand,  for  a  bipolar  outflow  the  absorption  velocities  may  be 
much  more  restricted  to  the  radial  ejection  speed  along  the  line  of  sight.  Conse¬ 
quently,  until  luture  observations  reveal  details  of  the  geometry  of  tile  nova  remnant, 
we  will  simply  assume  that  the  measured  radial  velocity  is  the  same  as  the  radial 
outflow  speed  in  our  direction  and  that  this  is  also  the  same  as  the  transverse  expan¬ 
sion  speed,  as  observed  in  the  plane  of  the  sky. 

31.  ten  Brummelaar.T.  A.  ef  a/.  The  CLASSIC/CLIMB  Beam  Combiner  atthe  CHARA 
Array .J.  Astron.  Instrum.  2,  1340004  (2013). 

32.  Monnier.J.D.ela/.  Michigan  Infrared  Combiner  (MIRC):  commissioning  results  at 
the  CHARA  Array.  Proc.  SPIE  6268, 62681P  (2006). 

33.  Bonneau,  D.eta/.  SearchCal:  a  virtual  observatory  tool  for  searching  calibrators  in 
optical  long  baseline  interferometry.  I.  The  bright  object  case.  Astron.  Astrophys. 
456,  789  (2006). 

34.  Boyajian.T.S.  etal.  Stellar  diameters  and  temperatures.  III.  Ma  n-sequence  A,  F,  G, 
and  K  stars:  additional  high-precision  measurements  and  empirical  relations. 
Astrophys.  J.  771, 40  (2013). 

35.  Crepp.J.R.  eta/.  The  dynamical  mass  and  three-dimensional  orbit  of  HR7672B:  a 
benchmark  brown  dwarf  with  high  eccentricity.  Astrophys.  J.  751,  97  (2012). 

36.  Skrutskie.  M.  F.etal.  The  Two  Micron  All  Sky  Survey  (2MASS).  Astron.  J.  131, 
1163-1183  (2006). 

37.  Henden,  A.  A.  The  AAVSO  International  Database,  http://www.aavso.org  (2013). 

38.  T uthill,  P.  et  al.  Double-Fourier  spatio-spectral  decoding.  In  Advances  in  Stellar 
Interferometry,  SPIEConf.  Ser.  6268,  62680X  (2006),  http://dx.doi.org/10, 1117/ 
12.672277. 

39.  Markwardt,  C.  B.  in  Astronomical  Data  Analysis  Software  and  Systems  XVIII,  Astron. 
Soc.  Pacif.  Cont  Ser.  411,  251-254  (2009). 

40.  Press,  W.  H„  Teukolsky,  S.  A.,  Vetterling,  W.  T.  &  Flannery,  B,  P.  Numerical  Recipes 
in  C  (Cambridge  Univ.  Press,  1992), 

41.  Wilson,  J.C.  etal.  In  Ground-based  Instrumentation  lor  Astronomy,  SPIE  Coni.  Ser. 
5492,  1295-1305  (2004);  http://dx.doi.org/10.1117/12.550925. 

42.  Herter,  T.  L.  ef  al.  in  Ground-based  and  Airborne  Instrumentation  for  Astronomy  II,  SPIE 
Coni.  Ser.  7014,  70140X  (2008);  http://dx.doi.org/ 1 0.1 1 17/12.789660. 

43.  Cushing,  M.  C.  Vacca,  W.  D,  &  Rayner,  J,  T.  Spextool :  a  spectral  extraction  package 
for  SpeX, a  0.8-5.5  micron  cross-dispersed  spectrograph.  Proc.  Astron.  Soc.  Pacif. 
116, 362-376(2004). 

44.  Muirhead,  P.S.ef  a/.  Precise  stellar  radial  velocities  of  an  M  dwarf  with  a  Michelson 
interferometer  and  a  med;um-resolution  near-infrared  spectrograph,  Publ.  Astron. 
Soc.  Pacif.  123,  709-724  (201 1), 

45.  Vacca,  W.  D„  Cushing,  M.  C,  &  Rayner,  J.  T.  A  method  of  correcting  near-infrared 
spectra  for  telluric  absorption.  Publ.  Astron.  Soc.  Pacif.  115,  389-409  (2003). 

46.  Banerjee,  D.  P.  K.,  Ashok,  N,  M„  Josh:,  V.  &  Evans,  A  Ongoing  near-infrared 
observations  of  V339  Del  (Nova  Del  2013).  Astron.  Telegr.  5404, 1  (2013). 

47.  Das.  R.  K„  Banerjee,  D.  P.  K„  Ashok,  N.  M.  &  Chesneau,  0.  Near-infrared  studies  of 
VI 280  Sco(Nova  Scorpii  2007).Mon.  Not  R.  Astron.  Soc.  391, 1874-1886  (2008). 

48.  Hjellming,  R.  M.  Wade,  C.  M„  Vandenberg,  N.  R.  &  Newell,  R.T,  Radio  emission  from 
nova  shells.  Astron  J.  84,  1619-1631  (1979). 


©2014  Macmillan  Publishers  Limited.  All  rights  reserved 


RESEARCH 


0  10  20  30  40  50 

Time  since  outburst  (days) 


Extended  Data  Figure  1  j  Integrated  emission  line  to  continuum  flux  ratios  emission.  The  down  turn  in  the  curve  during  die  last  measurement  is 

measured  from  infrared  spectroscopy.  The  squares  represent  the  H  band  probably  caused  by  a  rising  contribution  of  the  continuum  owing  to  the 

ratios  while  the  diamonds  represent  the  K'-band  ratios.  The  rise  in  the  emission  formation  of  dust, 

line  flux  is  consistent  with  an  increasing  contribution  from  optically  thin 
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Extended  Data  Figure  2  |  Closure  phases  measured  with  MIRC  on  ut  2013 
August  21.  The  measurements  are  plotted  against  the  length  of  the  maximum 
baseline  (B)  for  the  group  of  three  telescopes  that  were  used  to  form  the 
closure  phase.  The  non  zero  closure  phases  indicate  an  asymmetry  in  the  light 


distribution  that  is  not  point  symmetric.  The  colours  of  the  symbols  are  used  to 
differentiate  the  measurements  from  each  grouping  of  three  telescopes  in  the 
array.  The  closure  phases  are  measured  in  eight  wavelength  channels.  Error 
bars  represent  1<7  measurement  uncertainties. 
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Extended  Data  Figure  3  |  Example  of  M1RC  observations  obtained  on  ut  measurements,  colour-coded  to  match  the  baselines  on  the  left.  The  solid 
2013  August  21.  The  light  is  dispersed  over  eight  wavelength  channels  in  the  line  shows  the  best-fitting  uniform  disk  model.  The  small  black  dots  show  the 

H-  band,  a,  Coverage  of  the  interferometric  baselines  projected  on  the  plane  of  best  fitting  uniform  ellipse  modeL  Error  bars  represent  i  a  measurement 

the  sky  in  right  ascension  (RA)  and  declination  (Dec.)  in  units  of  spatial  uncertainties, 

frequencies  (u  =  BJX,  v=  B,/X).  b,  Squared,  normalized  visibility  amplitude 
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Extended  Data  Figure  4  |  Interferometric  visibilities  of  Nova  Del  201 3 
measured  with  the  CHARA  Array.  The  red  line  shows  the  best- fitting  model 
for  a  uniformly  bright,  circular  disk.  The  time  since  the  explosion  (in  dare)  is 


indicated  in  the  upper  right  comer  of  each  panel.  The  measurements  were 
obtained  with  the  CLASSIC,  CLIMB,  and  M1RC  beam  combiners  (see 
Extended  Table  1 ).  Error  bars  represent  la  measurement  uncertainties. 
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Extended  Data  Figure  5  |  Time  evolution  of  the  two-component  model  of  our  measurement  that  the  ring  contributes  an  average  of  68%  of  the  light  on  the 

Nova  Del  2013.  The  model  consists  of  a  circular  core  surrounded  by  a  halo  last  five  nights  (dust  emission  in  the  outer  layers).  Each  panel  is  12  mas  on  a 

ring.  The  expansion  rate  and  size  ratio  were  determined  by  minimizing  /'  side.  We  scaled  the  model  flux  by  the  infrared  magnitude  measured  on  each 

across  all  of  the  nights  while  allowing  the  flux  ratio  to  vary  night  by  night.  night  to  show  how  the  surface  brightness  changes.  The  time  since  the  explosion 

In  plotting  the  images,  we  used  the  flux  ratio  measured  directly  for  the  first  three  (in  days)  is  indicated  in  each  panel.  Intensity  refers  to  the  flux  per  unit  area, 

nights,  the  linear  fit  plotted  in  Fig.  3  of  the  main  paper  for  (  =  4-27  days,  and 
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Extended  Data  Table  1  |  Journal  of  observations  and  angular  diameters  measurements  of  Nova  Del  2013 


UT  Date 

MJD 

t  (days) 

Comb* 

Configuration 

Filter1 

Cal* 

OuD  (mas) 

1 

Xt 

N(l/2) 

2013  Aug  15 

56519.231 

1.0 

CL 

S1E2W1 

H,  K' 

1,2 

0.414  ± 

0.076 

1.19 

18 

2013  Aug  16 

56520.303 

2.0 

CL 

S1E1 

H 

1,2 

0.537  ± 

0.028 

0.23 

9 

2013  Aug  17 

56521.260 

3.0 

M 

W1S2S1E1E2W2 

H 

3 

0.652  ± 

0.005 

1.40 

180 

2013  Aug  18 

56522.284 

4.0 

C 

S1E1 

H 

1 

0.777  ± 

0.038 

1.77 

4 

2013  Aug  19 

56523.248 

5.0 

M 

W1S2S1E1E2W2 

H 

4,3 

0,863  ± 

0.004 

4.10 

180 

2013  Aug  20 

56524.325 

6.0 

CL 

E2W2W1 

H 

1 

1.096  ± 

0.056 

0.82 

6 

2013  Aug  21 

56525.220 

6.9 

M 

W1S2S1E1E2W2 

H 

4,3 

1.250  ± 

0.006 

11.5 

240 

2013  Aug  22 

56526.216 

7.9 

C 

E1W1 

K' 

4 

1.262  + 

0.046 

0.17 

2 

2013  Aug  23 

56527.310 

9.0 

CL 

S1E1W1,  S2E2W2 

H,  K' 

1,4,5 

1.643  ± 

0.077 

5.10 

15 

2013  Aug  24** 

56528.242 

10.0 

M 

W1S2S1E1E2W2 

H 

3 

1.885  ± 

0.016 

4.86 

32 

2013  Aug  26'* 

56530.208 

11.9 

M 

W1S2S1E1E2W2 

H 

3 

2.314  ± 

0.037 

1.09 

61 

2013  Aug  28 

56532.191 

13.9 

CL 

W1W2 

H,  K' 

4 

2.62  ± 

0.27 

2.49 

3 

2013  Aug  29 

56533.235 

15.0 

CL 

W1W2 

K' 

4 

3.133  ± 

0.076 

0.29 

2 

2013  Aug  30 

56534.273 

16.0 

CL 

E1E2,  S1S2 

K' 

1,4 

3.89  ± 

0.72 

0.19 

3 

2013  Sep  02 

56537.146 

18.9 

C 

E1E2 

K' 

4 

3.51  ± 

0.31 

1.10 

2 

2013  Sep  04 

56539.178 

20.9 

c 

E1E2 

K' 

1,4 

4.22  ± 

0.22 

3.32 

3 

2013  Sep  05 

56540.164 

21.9 

c 

E1E2 

K' 

1,4 

4.57  ± 

0.18 

0.04 

3 

2013  Sep  06 

56541.280 

23.0 

c 

E1E2 

K' 

1,4 

4.62  ± 

0.13 

0.40 

16 

2013  Sep  07 

56542.158 

23.9 

c 

E1E2 

K' 

1 

4.94  ± 

0.11 

0.23 

6 

2013  Sep  08 

56543.150 

24.9 

CL 

E1E2 

K' 

1 

5.11  ± 

0.12 

1.68 

5 

2013  Sep  09 

56544.156 

25.9 

CL 

E1E2 

K' 

1 

5.49  ± 

0.11 

1.76 

6 

2013  Sep  10 

56545.139 

26.9 

CL 

E1E2 

K' 

1 

6.09  ± 

0.17 

0.38 

2 

2013  Sep  13 

56548.160 

29.9 

CL 

S1S2 

K' 

1 

5.19  ± 

0.97 

0.31 

3 

2013  Sep  16 

56551.137 

32.9 

C 

S1S2 

K' 

1 

9.12  ± 

0.74 

0.85 

3 

2013  Sep  19 

56554.138 

35.9 

C 

S1S2 

K’ 

1 

9.33  ± 

0.30 

0.12 

5 

2013  Sep  23 

56558.121 

39.8 

C 

S1S2 

K' 

1 

10.96  ± 

0.18 

0.52 

5 

2013  Sep  26 

56561.159 

42.9 

C 

S1S2 

K' 

1 

11.83  ± 

0.98 

0.01 

2 

The  angular  diameters  were  computed  by  fitting  a  circular  uniform  disk  to  the  visibilities  from  each  night 
*Beam  combiner  codes:  C- CLASSIC,  CL  =  CLIMB,  M  =  MIRC. 

+  CLASSIC/CLIMB  H  filter:  /.  =  1.6731  pm,  1/  =  02854  pm:  CLASSIC/CLIMB  K'  filter: ;  =  2.1329  pm.  ix  =0.3489  pm;  MIRC  H  filter:  a  =  1.5-1 .7  pm  with  \/  =-0.034 pm  for  each  of  the  eight  spectral  channels. 
X  Calibrator  codes:  1  »  HD  190993  ((?  =  0.161  ±0.011  mas;  ref. 33),  2  =  HD  1 96740  (0  =  0.1 80  ±  0.01 2  mas;  ref.  33),  3  *  HD  206860(0  =  0.530  -0.0 15  mas;  ref.  34),  4  =  HD  190406(0=  0.584  ±  0.010  mas; 
ref.  35),  5  -  HD  195810  (&- 0.301  -0.021  mas,  ref.  33). 

§0nly  E2W1W2  arid  SI S2  were  recorded  on  ut  2013  August  24. Only  W1W2,  E1E2,  SI $2, and  E2W2  were  recorded  on  ut  2013  August  26. 
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Extended  Data  Table  2  |  Uniform  ellipse  models  of  Nova  Del  2013 


UT  Date 

MJD 

t  (days) 

Comb 

^major  (f"T13s) 

^minor  (m3s) 

PA  (c) 

\t 

2013  Aug  16 

56520.303 

2.0 

CL 

0.631  +  0.282 

0.503  ±  0.243 

81.8  ±  31.9 

0.14 

2013  Aug  17 

56521.260 

3.0 

M 

0.676  ±  0.011 

0.633  ±  0.007 

114.8  +  10.6 

1.30 

2013  Aug  19 

56523.248 

5.0 

M 

0.936  ±  0.005 

0.826  ±  0.004 

148.7  ±  2.1 

1.53 

2013  Aug  21 

56525.220 

6.9 

M 

1.291  +  0.007 

1.198  ±  0.005 

157.4  ±  3.5 

6.95 

2013  Aug  23 

56527.310 

9.0 

CL 

1.716  +  0.094 

1.535  ±  0.044 

138.5  ±  16.1 

2.49 
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Extended  Data  Table  3  |  Infrared  magnitudes  of  Nova  Del  2013 


MJD 

Fill 

Magnitude 

MJD 

Filt 

Magnitude 

MJD 

Filt 

Magnitude 

56519.179 

H 

5.12 

dr 

0.16 

56532.170 

K' 

4.44 

± 

0.23 

56542.149 

K' 

5.37 

± 

0.16 

56519.193 

H 

5.19 

± 

0.17 

56532.191 

K' 

4.53 

± 

0.16 

56542.158 

K' 

5.37 

± 

0.16 

56519.216 

H 

5.08 

± 

0.15 

56532.214 

H 

4.65 

± 

0.38 

56542.167 

K' 

5.35 

± 

0.15 

56519.231 

H 

5.20 

± 

0.14 

56533.223 

K' 

4.65 

± 

0.21 

56542.176 

K' 

5.32 

± 

0.16 

56519.250 

K' 

4.97 

± 

0.17 

56533.235 

K' 

4.60 

± 

0.20 

56543.125 

K' 

5.42 

i 

0.22 

56519.261 

K' 

5.04 

± 

0.17 

56534.249 

K' 

4.74 

i 

0.16 

56543.135 

K' 

5.44 

dr 

0.20 

56519.273 

K' 

4.93 

± 

0.16 

56534.265 

K' 

4.78 

± 

0.17 

56543.150 

K' 

5.42 

± 

0.19 

56520.187 

H 

4.20 

± 

0.20 

56534.273 

K' 

4.71 

± 

0.16 

56543.162 

K' 

5.36 

± 

0.19 

56520.218 

H 

4.24 

± 

0.19 

56534.404 

K' 

4.49 

± 

0.27 

56543.175 

K' 

5.30 

± 

0.18 

56520.247 

H 

4.16 

± 

0.17 

56537.130 

K' 

5.08 

± 

0.17 

56544.123 

K' 

5.32 

± 

0.22 

56520.278 

H 

4.21 

± 

0.16 

56537.146 

K' 

4.79 

-L 

0.18 

56544.132 

K' 

5.44 

± 

0.19 

56520.303 

H 

4.15 

± 

0.15 

56537.157 

K' 

4.46 

± 

0.31 

56544.145 

K' 

5.52 

± 

0.19 

56520.349 

H 

4.08 

± 

0.16 

56539.168 

K' 

5.10 

0.19 

56544.156 

K' 

5.52 

± 

0.17 

56520.375 

H 

4.11 

± 

0.16 

56539.179 

K' 

5.19 

it 

0.22 

56544.170 

K' 

5.42 

± 

0.15 

56520.412 

H 

4.07 

± 

0.18 

56539.190 

K’ 

5.62 

± 

0.23 

56544.184 

K' 

5.37 

± 

0.17 

56520.432 

H 

4.04 

± 

0.19 

56540.148 

K' 

5.25 

rh 

0.20 

56545.128 

K' 

5.25 

± 

0.29 

56521.241 

H 

4.09 

± 

0.12 

56540.164 

K' 

5.22 

± 

0.18 

56545.139 

K' 

5.46 

± 

0.30 

56521.263 

H 

4.05 

± 

0.12 

56540.179 

K' 

5.17 

± 

0.17 

56545.153 

K> 

6.56 

± 

0.54 

56522.265 

H 

3.74 

± 

0.15 

56540.194 

K' 

5.14 

± 

0.18 

56548.144 

K' 

5.38 

dr 

0.21 

56522.275 

H 

3.82 

0.14 

56540.212 

K' 

5.18 

± 

0.16 

56548.160 

K' 

5.26 

d= 

0.21 

56522.284 

H 

3.75 

± 

0.14 

56541.143 

K' 

5.29 

± 

0.15 

56548.171 

K' 

5.09 

dr 

0.22 

56522.295 

H 

3.70 

± 

0.14 

56541.160 

K' 

5.31 

± 

0.15 

56551.119 

K' 

5.43 

± 

0.27 

56523.219 

H 

3.63 

± 

0.15 

56541.178 

K' 

5.34 

± 

0.15 

56551.137 

K' 

5.49 

± 

0.26 

56523.248 

H 

3.46 

± 

0.16 

56541.195 

K' 

5.31 

± 

0.14 

56551.154 

K' 

5.28 

dr 

0.22 

56524.295 

H 

3.46 

± 

0.16 

56541.213 

K' 

5.22 

± 

0.14 

56554.118 

K' 

5.00 

± 

0.25 

56524.325 

H 

3.52 

± 

0.23 

56541.231 

K' 

5.30 

± 

0.12 

56554.129 

K' 

5.08 

± 

0.23 

56525.189 

H 

3.79 

± 

0.13 

56541.247 

K' 

5.29 

± 

0.13 

56554.139 

K' 

5.08 

± 

0.22 

56525.220 

H 

3.80 

± 

0.14 

56541.265 

K' 

5.34 

± 

0.13 

56554.150 

K' 

5.11 

± 

0.22 

56526.207 

K' 

3.58 

± 

0.15 

56541.281 

K' 

5.30 

± 

0.12 

56554.159 

K' 

5.17 

± 

0.22 

56526.224 

K' 

3.64 

± 

0.13 

56541.299 

K' 

5.22 

± 

0.14 

56558.102 

K' 

4.77 

± 

0.22 

56527.258 

K' 

3.86 

± 

0.13 

56541.315 

K' 

5.10 

± 

0.14 

56558.112 

K' 

4.80 

dr 

0.20 

56527.283 

K' 

3.68 

± 

0.13 

56541.333 

K' 

5.41 

± 

0.12 

56558.121 

K' 

4.77 

dr 

0.20 

56527.310 

K' 

3.84 

± 

0.16 

56541.350 

K' 

5.19 

± 

0.13 

56558.131 

K' 

4.79 

± 

0.20 

56527.324 

K' 

3.67 

± 

0.16 

56541.367 

K' 

4.86 

± 

0.15 

56558.139 

K' 

4.81 

± 

0.19 

56527.394 

H 

4.29 

dr 

0.15 

56541.389  ' 

K' 

5.20 

± 

0.14 

56561.129 

K' 

4.87 

± 

0.35 

56528.241 

H 

3.91 

± 

0.16 

56541.409 

K' 

4.49 

it 

0.17 

56561.159 

K' 

4.54 

± 

0.28 

56530.203 

H 

4.35 

± 

0.12 

56542.126 

K' 

5.38 

± 

0.22 

56530.263 

H 

4.54 

i 

0.11 

56542.141 

K' 

5.29 

± 

0.17 
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Extended  Data  Table  4  [  Effective  bandwidth  oftheCLASSIC/CLIIVlB 
observations 


UT  Date 

MJD 

t  (days) 

Bandwrdth(Nova) 

Bandwidth(Cal) 

2013  Aug  15 

56519.231 

1.0 

0.916 

±  0.295 

2013  Aug  16 

56520.303 

2.0 

1.015 

±  0.256 

2013  Aug  22 

56526.216 

7.9 

0.970 

±  0.193 

2013  Aug  28 

56532.191 

13.9 

1.023 

±  0.134 

2013  Aug  29 

56533.235 

15.0 

0.905 

±  0.107 

2013  Aug  30 

56534.273 

16.0 

0.879 

±  0.151 

2013  Sep  02 

56537.146 

18.9 

0.863 

±  0.104 

2013  Sep  04 

56539.178 

20.9 

0.844 

±  0.041 

2013  Sep  05 

56540.164 

21.9 

0.844 

±  0.121 

2013  Sep  06 

56541.280 

23.0 

0.876 

±  0.208 

2013  Sep  07 

56542.158 

23.9 

0.879 

±  0.117 

2013  Sep  08 

56543.150 

24.9 

0.817 

±  0.078 

2013  Sep  09 

56544.156 

25.9 

0.878 

±  0.086 

2013  Sep  10 

56545.139 

26.9 

0.792 

±  0.166 

2013  Sep  13 

56548.160 

29.9 

0.836 

±  0.098 

2013  Sep  16 

56551.137 

32.9 

0.933 

±  0.194 

2013  Sep  19 

56554.138 

35.9 

0.883 

±  0.162 

2013  Sep  23 

56558.121 

39.8 

0.902 

±  0.111 

2013  Sep  26 

56561.159 

42.9 

1.165 

±  0.492 

We  present  the  ratio  of  the  effective  bandwidth  of  the  interferometric  observations  of  Nova  Del2013 
compared  with  the  calibrator  stars. 
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